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摘 要 :采用 几何 精确 模型 与 阶 次 截断 方法 描述 叶片 的 空间 变形 ,建立 了 考虑 科 氏 效应 与 离心 刚 化 

的 叶片 - 桨 载 - 轴 耦 合 系统 动力 学 模型 ,推导 了 耦合 系统 振动 微分 方程 。 采 用 Rayleigh-Ritz 法 进 

寸 行 数值 离散 ,计算 了 不 同 刚度 比 、 半 径 比 转速 和 叶片 数 下 , 科 氏 效应 对 系统 模 态 频率 的 影响 。 研 究 
:叶片 - 桨 载 - 轴 耦 合 系统 模 态 中 包含 耦合 和 非 耦 合 两 类 模 态 ,耦合 模 态 频率 在 考虑 科 氏 效应 后 

低 ; 随 着 半径 比 增加 , 科 氏 效应 对 第 1 阶 耦 合 模 态 频 率 的 影响 逐渐 减弱 ,对 第 2 
殴 态 频率 的 影响 逐渐 增强 ; 转速 的 增加 会 提升 耦合 模 态 对 科 氏 效应 的 敏感 性 ;刚度 比 处 于 10 ~ 

人 Re eh ne 模 态 振 型 转换 后 科 氏 效应 的 影响 更 为 显著 ; 随 着 叶片 ww 

的 增加 , 科 氏 效应 对 模 态 频率 的 影响 逐渐 减弱 。 研 究 结 果 可 为 旋转 叶片 - 奖 载 - 轴 耦 合 系统 的 动力 学 

设计 提供 理论 参考 。 
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Coriolis effect on vibration characteristics of 
blade-hub-shaft coupling system 
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Abstract: Geometry exact model and order reservation approach were adopted to describe the spatial de- 
formation of the blade ,the dynamic model of the blade-hub-shaft coupling system considering the Coriolis 
effect and the centrifugal stiffening effect was established, and the vibration differential equation of the 
coupling system was derived. The Rayleigh-Ritz method was used for numerical dispersion ,and the influ- 
ence of the Coriolis effect on the system modal frequency was calculated under different stiffness ratios ,ra- 
dius ratios ,rotation speeds and blade numbers. The research shows that the modes of blade-hub-shaft cou- 
pling system consist of two types of coupled and uncoupled modes. The frequency of the coupled mode is 


significantly reduced after considering the Coriolis effect. As the radius ratio increases ,the influence of the 
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Coriolis effect on the first order coupled modal frequency is gradually weakened ,and the influence on the 


second-order coupled modal frequency is gradually strengthened; the raise of the rotational speed will im- 


prove the sensitivity of the coupled modal to the Coriolis effect; the blade modes shape conversion phe- 


nomenon occurs when the stiffness ratio is within the range of 10 -10” ,the effect of the Coriolis effect be- 


comes more significant after the mode shape conversion. As the number of blades increases ,the influence 


of the Coriolis effect on the modal frequency gradually decreases. The research results provide theoretical 


reference for the dynamic design of the rotating blade-hub-shaft coupling system. 


Key words :blade-hub-shaft coupling system; natural frequency; geometry exact; Coriolis effect; Ray- 


leigh-Ritz method 


旋转 叶片 - 桨 载 - 轴 结 构 常见 于 各 种 动力 与 推进 
设备 ,叶片 在 旋转 时 工 况 复杂 ,其 振动 特性 一 直 备 受 
关注 。 国 内 外 研究 者 对 旋转 叶片 结构 进行 了 大 
量 研究 5 。 
亏 刚 柔 耦 合 模型 是 一 种 应 用 较为 广泛 的 旋转 粱 动 
为 学 特性 分 析 模型 。 杨 辉 等 ”从 刚 柔 耦 合 角度 针对 
旋 团 柔性 梁 - 刚 性 转轴 系统 进行 研究 ,利用 一 次 近似 
模型 对 旋转 梁 的 动力 学 响应 并 进行 了 仿真 计算 及 实 
骏 验 证 ,表明 旋转 粱 需 考虑 离心 刚 化 作用 。 吴 胜 宝 
等 中 在 计算 旋转 梁 动 力学 特性 时 在 纵向 变形 位 移 中 
证 及 了 横向 弯曲 引起 的 轴 向 缩短 ,表明 柔性 梁 中 既 
可 得 在 动力 刚 化 效应 ,也 可 存在 动力 柔 化 效应 。 郭 
4 涡 等 ”采用 哈密 尔 顿 原理 建立 了 含 偶 应 力 弹性 体 
作 定 轴 刚 体 转动 的 运动 与 变形 耦合 动力 学 模型 , 针 
对 天 同 的 角速度 变化 规律 ,讨论 了 旋转 悬臂 方 梁 的 
非 线 性 动力 学 响应 特征 。 盛 国 刚 等 “建立 了 支承 运 
动情 况 下 旋转 粱 的 刚 - 柔 耦 合 振动 非 线 性 动力 学 模 
型 E 分 析 了 刚性 角 位 移 与 梁 的 动力 学 响应 之 间 的 关 
系 GDAn 等 ”分 析 了 变 刚 度 的 旋转 梁 瞬 态 动力 学 响 
应 问题 。 文 献 [8-9] 利 用 梁 中 心 线 的 倾斜 角 和 拉 伸 
应 变 来 描述 欧 拉 伯 努 利 梁 的 变形 ,分 析 了 具有 恒定 
角速度 情况 下 的 固有 频率 及 振 型 ,表明 在 考虑 弯曲 
和 拉 伸 耦合 效应 时 , 会 出 现 频率 曲线 转向 现象 。 
Zhao 等 '" 利用 牛顿 定律 推导 出 梁 的 弹性 运动 和 兹 
载 转动 的 控制 方程 ,分 析 了 旋转 欧 拉 - 伯 努 利 梁 的 稳 
态 振动 ,讨论 了 转轴 转动 惯量 与 半径 对 旋转 梁 频 率 
与 振 型 的 影响 。 刚 柔 耦 合 模型 能 够 分 析 柔 性 粱 在 特 
定 转动 规律 下 的 动力 学 特性 ,但 实际 中 的 转轴 并 不 
能 够 完全 等 效 为 刚体 。 为 了 更 加 精确 地 分 析 叶 片 - 
桨 载 -转轴 系统 动力 学 特性 , 需要 考虑 转轴 的 弹 
性 振动 。 

一 些 学 者 考虑 了 转轴 的 弹性 变形 ,分 析 了 叶片 - 
转轴 耦 合 系统 动力 学 特性 。 贺 天 月 等 … 针对 共 轴 


直升机 旋 波 系统 进行 了 多 体 动力 学 分 析 , 模 型 中 包 
含 了 柔性 旋 波 与 柔性 轴 , 利 用 多 体 动力 学 软件 AD- 
AMS 进行 了 仿真 计算 并 进行 了 实验 验证 。 结 果 表 
明 , 考 虑 了 旋翼 与 轴 柔 性 的 计算 结果 与 实验 数据 吻 
合 良好 。 徐 自力 等 ”利用 对 多 级 叶 盘 - 轴 耦 合 结构 
的 振动 特性 进行 了 计算 分 析 , 发 现 考虑 轴 弹 性 扭转 
变形 作用 后 ,系统 模 态 中 增加 了 新 的 耦合 模 态 ,耦合 
模 态 振动 特征 与 单纯 叶片 振动 特征 存在 明显 区 别 。 
朱 由 锋 等 ”基于 有 限 差分 法 建立 了 水 平 旋转 梁 的 
动力 学 模型 ,分 析 了 不 同 转速 下 的 弹性 梁 固 有 频率 
变化 规律 ,但 未 考虑 系统 旋转 带 来 的 科 氏 效应 。 李 
朝 峰 等 “建立 了 包含 多 个 弹性 叶片 - 轮 盘 -转轴 耦合 
系统 ,发现 了 叶片 数 、 转 盘 位 置 对 耦合 系统 固有 频率 
均 存 在 明显 影响 ,但 并 未 深入 分 析 科 氏 效 应 对 系统 
动力 学 特性 的 影响 。 

对 于 旋转 叶片 , 科 氏 效应 对 其 振动 特性 具有 显 
著 影 响 ,尤其 当 叶 片 柔性 较 大 时 ,叶片 变形 非 线性 程 
度 增 加 。 为 了 更 为 精确 地 分 析 科 氏 效 应 对 叶片 - 浆 
坑 - 轴 耦合 系统 动力 学 特性 的 影响 , 需 准确 描述 叶片 
变形 。 本 研究 利用 几何 精确 梁 模 型 并 结合 阶 次 截断 
方法 ,考虑 了 叶片 拉 伸 、 弯 曲 及 扭转 的 三 维 空间 运 
动 ,建立 了 考虑 科 氏 效应 与 离心 刚 化 作用 的 高 精度 
叶片 - 桨 坑 - 轴 耦 合 动力 学 模型 ,并 通过 数值 模拟 研 
究 了 不 同 半径 比 刚 度 比 . 叶 片 数目 下 科 氏 效应 对 耦 
合 系统 振动 特性 的 影响 规律 。 


1 叶片 - 桨 席 - 轴 耦 合 动力 学 模型 


叶片 - 桨 载 - 轴 耦 合 系统 包含 柔性 叶片 .刚性 桨 
坑 以 及 可 发 生 弹 性 扭转 变形 的 转轴 。 如 图 1 所 示 ， 
OXoJ2 为 惯性 参考 系 ,坐标 轴 指 向 遵循 右手 定 则 ， 
0 在 系统 转轴 中 心 线 上 ,2 轴 垂 直 指向 面 外 。 整 个 
系统 以 角速度 w 绕 Z 轴 旋 转 。 定 义 桨 载 随 体 旋转 


坐标 系 为 0X,7Z,Xi 与 Xo 夹 角 为 wi,1 为 转动 时 
间 。 轴 的 弹性 扭转 变形 使 得 浆 载 在 基础 的 旋转 角度 
上 附加 了 一 个 扭转 角 pg。 叶 片 随 体 坐 标 系 OX,Y,Z， 
% 与 叶片 未 变形 时 的 中 心 线 重合 ,RR 为 桨 融 半 径 。 


图 1 叶片 - 浆 载 - 轴 结 构 坐 标 系 统 描述 


Fig.1 Description of coordinate systems 


™ 
> of blade-hub-shaft structure 


村 图 2 所 示 为 叶片 随 体 坐标 系 中 变形 示意 图 。 图 

中 .Le 为 固定 在 叶片 截面 上 的 坐标 系 ;坐标 原点 

于 截面 中 心 ; 轴 与 叶片 中 心 线 重合 ,垂直 于 叶 

片 组 向 ;7 轴 平 行 于 叶片 弦 向 。 叶 片 变形 后 ,0, 沿 
E60 方向 的 位 移 分 别 为 .vww。 
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图 2 叶片 变形 示意 图 
Fig.2 Blade deformation diagram 

根据 上 述 坐 标 描述 ,叶片 未 变形 时 任意 位 置 处 

的 质点 可 表示 为 
Ro=[R+x 0 0] +[0 7 2] (1) 

式 中 x 为 奖 载 至 该 质点 在 郊 轴 上 的 投影 距离 。 

叶片 变形 后 任意 位 置 处 的 质点 可 表示 为 

Ri=[R+ix+tu v 2 +A.:[0 7 0 (2) 
式 中 4 为 欧 拉 变换 矩阵 。 为 了 简化 计算 ,对 矩阵 中 
各 元 素 进 行 阶 次 截断 ,保留 至 二 阶 ,得 到 简化 后 的 欧 


1 2 Ww 
4=| -0w-v 1-0wv 0 (3) 
—-w -0 -0-wwvw 1 
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式 中 :9 为 叶片 截面 转角 ; * ' 为 * 对 x 求 偏 导 。 
由 格林 应 变 准 则 ,根据 变形 前 后 的 位 置 矢量 可 
得 到 叶片 应 变 张 量 为 


1 
£1= 7 (GG -8 8;) (4) 
其 中 
0C， 9Go 证 
> st 三 5 
i [x 7 9 (5) 


忽略 二 阶 以 上 的 高 阶 项 , 式 (5) 中 应 变 张 量 可 
表示 为 


=C,V, (6) 
其 中 
1 T+ ¢ -7 
C, -| - 0 0 
0 刀 0 0 
2 + +w’)? 
V, = (0' +w'w)’ 
vO-w 
vw +w"0 
式 中 ,下 标 p 表示 势能 相关 项 。 
式 (6) 变 分 可 得 
Se =C,H.,50, (7) 
其 中 
Q,=[0 wu vv w 0 vw w”]"， 
「 0 0 vw] 
1 0 0 0 
v" 0 0 0 
H,=|jw 2(0 +w’w")v 0 0 
0 2(0 TWO 0 0 
0 2(0 +w'v)w 0 1 
LO 0 -1 0 


式 中 :@, 为 叶片 应 变 广义 坐标 向 量 ; 五 , 为 应 变 变 分 
转换 矩阵 。 

叶片 弹性 势能 变 分 由 叶片 应 变 及 其 变 分 项 
表示 为 

80, = se'Eedv + | | GpspySw "Edéds (8) 


式 中 : 脚 标 b 表示 叶片 相关 项 ;第 一 积分 项 表示 由 叶 

片 变形 产生 的 应 变 势 能 ;第 二 积分 项 为 叶片 旋转 时 

离心 力 造成 的 叶片 拉 伸 变形 势能 ;E 为 叶片 的 弹性 

模 量 和 矩阵 ;V 为 叶片 积分 微 元 体积 ;p, 为 叶片 密度 ;S 

为 叶片 截面 面积 ;w 为 系统 转速 ;mm 为 叶片 长 度 。 
将 式 (6) 和 式 (7) 代 入 式 (8) 中 ,整理 可 得 
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Lb 
8U, = | ,5QW HY S, Vode + 


Lb Lb 六 
| 5601 | PSw2Ededx 


) 
二 


S = | GEC,s, 
V =[0 1+ wv w 0 0 0] 
根据 DAlembert 原理 ,叶片 动能 变 分 表达 式 为 
67, = - | 8R' p, RdV (10) 
式 中 :R, 为 惯性 坐标 系 OXo82 中 的 质点 位 置 拓 量 ; 
及, 为 位 置 矢量 对 时 间 的 二 阶 导数 。 
叶片 随 体 坐标 系 OX,Y,Z 中 的 位 置 矢量 R, 与 惯 
性 邓 中 的 位 置 天 量 R。 之 间 满 足 如 下 关系 。 
R,=A,R, (11) 
式 申 : 4， 为 惯性 系 与 非 惯性 系 之 间 的 矢量 旋 


cos(wf+D) sin(wt+p) 0 


4,=| -sin(witp) cos(of+p) 0| (12) 
CU 矢量 旋转 矩阵 中 ,包含 了 系统 旋转 参数 wt 以 及 
桨 融 扭 转角 9 两 部 分 。 

3 对 式 (11) 取 变 分 ,位 置 矢量 变 分 SR, 可 表示 为 


5R, =A,C.H.8Q, +A,C. VS9 (13) 


其 串 
-sin(wl+op) cos(wit+wp) 0 
和 | -ee -sin(owf+o) 中 
0 0 0 
Ce=[5a 7 XT 《xDxs]， 
[rst+%+u | 
v 
10 
vy Ow’ 
V. = 1 
0 
Ov’ —w’ 
-0 
L 1 J 


[1 0 0 0 0 0 
0 10 0 0 0 
0 0 1 0 0 0 
0 0 0 -ww =-1 0 
H=I0 0 0 vw -60 -0 
0 0 0 1 0 0 
0 .00 7 0 -1 
0 0 0 -1l —w’ 一 2 
LO0 0 0 0 0 0 J 


式 中 :Q, 为 叶片 动能 广义 坐标 向 量 ; V 为 运动 参数 
向 量 ; 及 , 为 运动 变 分 转换 矩阵 。 
式 (10) 中 的 有 R, 可 由 式 (11) 左 右 两 端 对 时 间 求 
二 阶 导数 得 到 
R=A,C(V+2QNV +QV) (14) 


其 中 
0 w 0 -w op 0 
2=| -wo 0 02=| -yp -w 0 
0 0 0 0 0 0 


将 式 (13) 和 式 (14) 代 入 式 (10) , 则 叶片 动能 变 
分 57, 写 为 


Lb . 六 
387，= - 上 8 OH' SV +20V, + QV)dx- 


Lb 六 。 入 
| sp V' SV +20V, + QV,) dx 
0 


(15) 
其 中 


S, = p, | CT Cds,S, = p, | CT ATA, Cids 
s 3 


式 (15) 中 ,20 VY, 项 即 为 科 氏 效应 项 。 
类 似 于 式 (10) , 奖 载 动能 变 分 57, 可 表示 为 
57, = -pp (16) 
式 中 五 为 桨 载 的 转动 惯量 。 
同样 可 求 得 转轴 的 动能 变 分 57. 为 
37 =- 三 sp Wdz (17) 
式 中 :yy(z,?) 为 轴 在 Z 方向 上 的 扭转 角 分 布 函 数 ;p、 
为 轴 密 度 ;/. 为 轴 截 面 极 惯性 矩 ;元 为 轴 长 度 。 
轴 扭 转变 形 的 弹性 势能 变 分 5U. 为 
8U = [wep (18) 
式 中 C 为 转轴 的 切 变 模 量 。 
至 此 已 推导 出 叶片 - 浆 载 - 轴 耦 合 系统 各 部 分 的 
动能 及 势能 表达 式 。 根 据 Hamilton 原理 ,叶片 - 浆 


载 - 轴 和 耦合 系统 的 运动 满足 如 下 关系 。 
六 sw + 8U. -5 -5 - 87,)dt = 0 (19) 


式 中 :ti 各 为 任意 时 刻 。 
采用 Rayleigh-Ritz 法 对 运动 变量 进行 解 耦 ,得 
u(x,t) = D(x)q,.(i) ,v(x,t) = D(x)g,(i) 
w(x,i) = D(x)q,(i) ,0(x,i) = B(x) qo(t) 
水 (zt =B(z)é(t) ,p(t) = DB(L.)é(i) 
(20) 
式 中 :B(x) 为 已 知 的 形 函 数 向 量 , 可 根据 计算 精度 
与 效率 的 需求 设置 不 同 阶 数 ;9 , (1) 为 对 应 各 个 叶 
片 运 动 变 量 的 广义 坐标 向 量 ;E(t) 为 对 应 轴 扭 转运 
动 的 广义 坐标 向 量 。 需 要 指出 的 是 , 桨 红 的 扭转 角 
度 -p 满足 g(t) =yw(L,,t)。 
整体 广义 坐标 向 量 8 可 写作 
cS Q=[g, q, gq, 9 |] (21) 
YJ 结合 式 (20 ) 与 式 (21) ,各 广义 坐标 向 量 可 


1 


(22) 
y=NQ,p=NiQ 


式 虽 :入 , ~ N, 矩阵 为 对 应 各 个 运动 坐标 向 量 的 形 函 
闫 绚 阵 。 
-将 式 (22) 各 项 代入 式 (19) 可 得 
QO(F, -FPF,+F.,-F.) =0 (23) 


人 =NQ,0.=N,0 


式 中 
. Lib 本 Dh Lb 
FE | N' HY S, Vdx + | Ni | prSw Edédx, 
Fi 可 一 Nh 9p， 
Is Ls . 
F.= | NGIw'd + | No 人 yadz， 
0 0 


Lb , . 
F, = | NiH' SS,(V, +20V + QV)dx- 


Lb p a 。 号 
| NI VI SV +20V. + QV,)dx 
0 


式 (23) 对 于 任意 88 均 成 立 ,可 得 系统 隐 式 运动 微 
分 方程 , 即 


=F,-F,+F,-F,=0 (24) 
系统 质量 矩阵 M 刚度 矩阵 玉 和 科 氏 效应 和 矩阵 
G 可 分 别 表示 为 
01 917 oar 
so 30,G 0 (25 ) 
系统 振动 微分 方程 可 写作 
MO+GO+KO=0 (26) 


需要 注意 的 是 ,在 能 量 泛 函 了, 的 第 二 项 中 考虑 
了 离心 力 引 起 的 形变 势能 ,因此 ,利用 总 体能 量 泛 函 
求 得 的 刚度 矩阵 K 中 已 包含 了 离心 刚 化 作用 。 

对 于 式 (26) 的 复 特 征 值 问题 ,采用 状态 空间 方 
法 对 系统 进行 求解 引入 状态 空间 变量 


X=[Q0]" (27) 
根据 式 (27) 将 式 (26 ) 改 写 为 
MX+KX=0 (28) 


其 中 
~ 10 MI- -M0 
ural | 0 刘 
则 式 (28) 的 特征 值 问题 转 化 为 
[AM +K|=0 (29) 
由 式 (29) 即 可 求 得 叶片 - 桨 载 - 轴 系 统 特征 频率 
与 模 态 振 型 。 


2 算 例 分析 


2.1 模型 验证 


叶片 - 桨 载 - 轴 看 合 系统 的 几何 物理 参数 如 表 1 
所 示 。 
表 1 叶片 - 桨 载 - 轴 耦 合 模型 参数 
Tab.1 Parameters of blade-hub-shaft coupling model 


参数 名 称 参数 值 
叶片 截面 积 S/m? 2.5x10-5 
叶片 密度 ph/(kg :mm 一 ) 2 700 
转轴 密度 ps/(kg mm ) 2 700 
转轴 长 度 L,/m 0.3 
叶片 长 度 LA/m 1 
叶片 弹性 模 量 /Pa 7.1x10" 
转轴 切 变 模 量 6/Pa 2.67 x10" 


为 使 分 析 更 具 一 般 性 , 取 如 下 形式 的 量 纲 为 一 
的 参数 进行 讨论 。 
CU 


_ R PS 二 PS __ 
L, ee El, EL, We ElL. 


式 中 :o 为 桨 载 与 叶片 半径 之 比 ;x 为 系统 量 纲 归 一 
化 转速 ;y 为 系统 量 纲 归 一 化 频率 jw 为 叶片 弯曲 -得 
扭转 刚度 比 ;为 叶片 截面 惯性 矩 。 

上 述 量 纲 归 一 化 参数 给 定 后 , 桨 载 半 色 让 . 桨 载 
的 转动 惯量 、 轴 的 玲 面 极 惯性 矩 均 可 求 得 。 


532 


首先 ,对 本 研究 算法 的 准确 性 进行 验证 。 综 合 
考虑 计算 效率 与 精度 ,本 研究 使 用 的 形 函 数 阶 数 为 
15 ,计算 时 消去 科 氏 效应 矩阵 G 以 忽略 科 氏 效应 。 
表 2 给 出 了 利用 本 研究 提出 的 数值 模型 计算 的 
0 =0,k=% 时 , 均 质 旋转 叶片 的 前 2 阶 量 纲 归 一 化 
频率 y 与 已 有 文献 计算 结果 "的 对 比 情况 。 表 3 
给 出 了 考虑 科 氏 效应 不 同 转速 和 半径 比 下 叶片 基础 
频率 与 文献 结果 ”对 比 结果 。 可 以 看 出 ,本 研究 数 
值 计算 结果 与 现 有 文献 数据 吻合 良好 ,具有 和 较 高 的 
计算 精度 。 

表 2 不 考虑 科 氏 效应 时 ,均匀 叶片 在 不 同 转 速 下 的 
前 2 阶 频率 (o =0, 人 =oo) 


Tab.2 First 2 step frequencies of a uniform blade under 


various rotating speeds without Coriolis effect(o =0, 从 =oo) 
4 型 


?1 ?72 
文献 [ 10 ] 本 研究 文献 [ 10 ] 本 研究 
3.5160 3.5160 22.035 22.036 
3.543 4 3.543 4 22.159 22. 161 
4.077 8 4.077 9 24.948 24.950 
5.059 6 5.059 8 32.113 32. 120 


表 3 考虑 科 氏 效应 时 ,不 同 转速 和 半径 比 下 
的 叶片 基础 频率 (=o%) 
Fundamental frequency of a uniform blade under various 


ting speeds and radius ratios with Coriolis effect(j, = oo) 


K 文献 [9 1 本 研究 偏差 /% 
4 2 5.049 5 5.051 0 0.029 7 
10 13.349 7 13.356 7 0.052 4 
50 41.241 1 41.274 5 0.081 0 
5 2 7.097 4 7.094 1 —0.046 5 
10 27.467 0 27.5067 0.1445 
50 73.473 6 73.563 0 0.121 7 


图 3 给 出 了 量 纲 归 一 化 转速 x 由 0 增加 至 100 
时 , 量 纲 归 一 化 频率 y 的 变化 规律 。 可 以 看 出 ,在 转 
速 较 小 时 ,两 种 曲线 并 未 出 现 明显 差异 , 科 氏 效应 并 
不 显著 ; 当 转 速 增加 后 , 科 氏 效应 加 剧 ,两 种 曲线 发 
生 分 离 。 在 图 中 标注 了 前 4 阶 模 态 频率 曲线 分 别 记 
为 C-1 F-1、.C-2 F-2。 可 以 看 出 , 科 氏 效应 对 C-l 与 
C-2 模仿 频率 曲线 影响 较为 显著 ,而 对 F-1、F2 模 态 
频率 曲线 影响 较为 微弱 。 


最 纲 归 一 化 频率 》 


--- 不 考虑 和 成 
一 考虑 科 氏 多 


0 20 40 60 80 100 
量 网 归 一 化 转速 天 


到 3 科 氏 效应 对 叶片 浆 坑 - 轴 系 统 频率 的 影响 


Fig.3 Influence of Coriolis effect on natural frequencies 
of the blade-hub-shaft system 

图 4 中 给 出 了 上 述 4 阶 叶 片 - 浆 载 - 轴 耦 合 模 态 
振 型 。 可 以 看 出 ,C-l 模 态 为 叶片 1 阶 弯曲 与 轴 的 
扭转 耦合 模 态 ;F-1 模 态 为 单纯 叶片 弯曲 振动 ,未 与 
轴 的 扭转 振动 发 生 耦 合 ; C-2 模 态 为 叶片 的 2 阶 弯 
曲 与 轴 的 扭转 耦合 模 态 , 且 与 C-1 模 态 振 型 相 比 , 叶 
片 的 弯曲 方向 相反 ;F-2 模 态 为 单纯 叶片 2 阶 弯曲 模 
态 ,未 与 轴 扭 转发 生 耦 合 

C-l 与 C-2 模 态 均 为 叶片 弦 向 弯曲 与 轴 扭 转 耦 
合 模 态 。 从 科 氏 效应 产生 的 根源 上 分 析 , 叶 片 截面 
3 个 正 交 方 向 的 运动 变量 量 中 , 轴 向 位 移 wx 和 弦 向 位 
移 " 与 旋转 角速度 方向 正 交 ,会 产生 科 氏 力 ;而 叶片 
挥舞 方向 位 移 w 与 旋转 角速度 重合 ,不 会 产生 科 氏 
力 。 因 此 C-1、C-2 耦合 模 态 受到 的 科 氏 效应 影响 最 


{aye-l ihy FP-1 
1 -了 (di [-2 


图 4 耦合 系统 前 4 阶 模 态 振 型 图 (x =0) 
Fig.4 First 4 modal shapes of coupling system(k =0) 
上 述 数 值 结果 说 明 ,在 较 高 转速 下 ,叶片 的 面 内 
弯曲 振动 必须 考虑 科 氏 效应 的 影响 。 


文中 采用 频率 相对 偏差 玉 参数 表征 科 氏 效应 
的 强 弱 。 
E,=Y x100% (30) 
yi 


式 中 :7; 为 第 i 阶 考虑 科 氏 效应 的 模 态 频率 ;y; 为 多 
i 阶 不 考虑 科 氏 效应 的 模 态 频率 。 


站 


2.2 b, 随 半径 比 的 变化 规律 


图 5 给 出 不 同 半径 比 情况 下 ,前 2 阶 耦 合 模 态 
中 叶片 量 纲 归 一 化 振 型 的 变化 趋势 。 图 中 横 轴 表 示 
叶片 轴 向 位 置 , 纵 轴 为 "方向 上 的 模 态 振 型 分 量 。 
由 图 5 可 以 看 出 :半径 比 由 0 增加 至 0.5 时 ,C-1 阶 
模 态 振 型 的 变形 程度 依次 降低 ; 随 着 半径 比 增加 ,C- 
2 阶 模 态 叶片 模 态 振 型 的 变形 程度 依次 升 高 。 


1.0 


0.8 
炬 
守 0.6 
了 
六 04 
外 
本 = 
-4 02 04 06 08 1.0 0 02 04 06 08 1.0 
叶片 轴 向 位 置 wl 叶片 轴 癌 位置 只 有 
所 二 2 模 坊 
SS (a) C-1 模 态 (b) C-2 模 态 
a 图 5 不 同 半径 比 下 ,前 2 阶 耘 合 模 态 叶片 
2 量 纲 归 一 化 振 型 (1 =10') 
CH. 5 Normalized first 2 coupling mode shapes of blade 
CN with various radius ratios (人 = 10” ) 


民国 6 给 出 了 半径 比 在 0 ~0.5 范围 变化 时 的 前 2 
阶 查 合 模 态 频率 相对 偏差。 


K=20 
-~--.-K=40 
区 K=60 
一 --K=80 
[3 Ar 本 
20 上 + A 
一 
0 ~ 要 
带 一 -一 ” 
全 
族 10F see 
聚 
本 
I | 
0 0 2 037 Q4. (3 0 0&1 162 43 1 人 4 05 
半径 比 o 半径 比 o 
(a) C-1 模 态 (b) C-2 模 态 


图 6 不 同 半径 比 下 ,前 2 阶 耦 合 模 态 频率 
的 相对 偏差 人 =10” ) 


Fig.6 Relative error of first 2 coupling mode frequencies 


with various radius ratios (人 = 10” ) 


由 图 6(a) 可 以 看 出 : 当 k=20 时 ,C-l 模 态 频率 


相对 偏差 E; 始终 保持 在 10% 以 下 ; 当 =40 时 ,局 


随 着 半径 比 增加 逐渐 降低 ,这 种 趋势 在 k = 60 和 
K=80 情 况 下 更 为 明显 。 由 图 6(b) 可 以 看 出 ,C-2 
模 态 频率 显示 出 与 C-1 模 态 相反 的 规律 ;不 同 的 转 
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速 下 ,C-2 模 态 频率 相对 偏差 随 着 半径 比 的 增加 均 
有 一 定 程度 的 上 升 ,在 高 转速 下 的 上 升幅 度 更 大 。 
产生 上 述 现 象 的 原因 为 , 科 氏 效应 强度 与 叶片 vo 方 
向 变形 程度 相关 。 结 合 图 5 的 叶片 变形 情况 ,半径 
比 增加 后 ,C-1 模 态 中 的 叶片 v 方 向 的 变形 程度 减 
弱 , 科 氏 效应 随 之 降低 ;而 C2 模 态 中 的 叶片 "方向 
的 变形 程度 随 半 径 比 增加 逐渐 增 大 , 科 氏 效应 随 之 
提高 ;转速 提高 后 , 科 氏 效应 整体 增强 ,相对 偏差 进 
一 步 增 大 。 

上 述 数 值 分 析 结果 表明 , 科 氏 效应 对 不 同 模 态 
的 影响 规律 存在 差异 ,在 工程 分 析 过 程 中 需 单独 考 
虑 各 阶 模 态 下 的 频率 变化 规 得 


2.3 已 随 刚度 比 的 变化 规律 


图 7 给 出 了 不 同 刚度 比 情况 下 ,前 2 阶 模 态 中 
叶片 量 纲 归 一 化 振 型 的 变化 趋势 。 由 图 7(a) 可 以 
看 出 :刚度 比 由 10 增加 至 10 时 ,C-l 模 态 叶片 振 型 
的 变形 程度 依次 增 大 ,其 中 :刚度 比 由 10° 增 加 至 10 
时 ,叶片 模 态 振 型 变化 较为 明显 ;而 刚度 比 由 10” 增 
加 至 10 时 ,叶片 模 态 振 型 变化 很 小 。 由 图 7(b) 可 
以 看 出 ,刚度 比 为 10' ~ 10 时 的 C2 模 态 叶片 振 型 
曲线 接近 为 叶片 的 1 阶 弯曲 ,刚度 比 为 10” ~ 10 时 
的 叶片 振 型 曲线 接近 叶片 的 2 阶 弯 曲 。 这 说 明 在 上 
述 两 个 范围 之 间 存 在 一 个 叶片 模 态 转换 现象 , 模 态 
转换 前 后 C-2 模 态 叶片 振动 形式 发 生 明 显 变 化 。 


1.0 1.2 
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训 
窗 0.8 
六 0.6 
苹 0.6 
吕 0.4 
并 0.4 
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0 02 04 06 08 1.0 0 02 04 06 08 1.0 
叶片 办 向 位 置 WA 叶片 轴 问 位 置 WA 
(a) C-1 模 态 (b) C-2 模 态 
图 7 不 同 刚度 比 下 ,前 2 阶 耦 合 模 态 
叶片 量 纲 归 一 化 振 型 (o =0. 1) 
Fig.7 Normalized first 2 coupling mode shapes of blade 


with various stiffness ratios (0 =0.1) 
图 8 给 出 了 前 2 阶 耦合 模 态 频率 相对 侦 差 随 刚 
度 比 的 变化 规律 。 由 图 8(a) 可 以 看 出 ,C-1 模 态 频 
率 相 对 侦 差 曲线 呈 先 降低 后 增加 的 趋势 。 结 合 图 7 
分 析 : 当 刚度 比 在 10' ~10 范围 内 变化 时 ,叶片 变形 
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程度 较 低 , 科 氏 效 应 较 弱 ,但 由 于 此 时 C-l 模 态 频率 
本 映 较 低 ;因此 相对 偏差 较 大 。 随 着 刚度 比 的 增加 ， 
C-1 模 态 频率 迅速 升 高 ,而 科 氏 效应 并 未 明显 增强 ， 
因此 相对 偏差 逐渐 降低 。 刚 度 比 在 10” ~10 范围 内 
变化 时 , 叶 卢 变形 程度 逐渐 增 大 , 科 氏 效应 随 之 增 
强 , 因 此 相对 偏差 逐渐 升 高 。 刚 度 比 在 10” ~ 10 范 
围 内 时 ,叶片 变形 程度 基本 不 变 , C-1 模 态 频率 也 基 
本 保持 不 变 , 因 此 相对 偏差 维持 在 一 个 定 值 。 此 时 ， 
由 于 轴 刚 度 较 大 ,基本 不 发 生 扭 转 振动 , 故 半径 比 的 
变化 对 相对 偏差 影响 不 大 。 从 图 8(b) 可 以 看 出 , 刚 
度 比 在 10 ~10 区 域内 相对 偏差 出 现 明显 的 上 升 。 
由 前 文 分 析 可 知 , 该 区 域内 发 生 了 叶片 模 态 转换 , 模 
态 转换 前 叶片 变形 相对 较 小 , 科 氏 效应 相对 较 弱 , 相 
对 偏差 较 低 ; 发 生 模 态 转换 后 ,叶片 变形 增 大 , 科 氏 
效应 增强 ,相对 偏差 增 大 。 


相对 偏差 E19% 


频 识 


刚度 比 x 


(b) C-2 模 态 


(a) C-1 模 态 


图 [83 不 同 刚度 比 下 前 2 阶 灯 合 模 态 频率 相对 偏差 (« =40) 
年 g 8 Relative error of first 2 coupling mode frequencies 


iE . , 
with various stiffness ratio (K =40) 


上 〇 结果 表明 ,过 大 和 过 小 的 转轴 刚度 均 会 引起 明 
显 的 科 氏 效应 。 由 于 科 氏 效应 会 降低 系统 频率 , 引 
起 失 稳 ,因此 需要 在 系统 设计 时 将 转轴 与 叶片 的 刚 
度 比 控制 在 10” 附近。 


2.4 ;随时 片 数 的 变化 规律 


图 9 给 出 了 叶片 数目 入 为 1 ~5 时 的 前 2 阶 耦 
合 模 态 叶 片 量 纲 归 一 化 振 型 。 由 图 9 可 以 看 出 , 随 
着 叶片 数目 增加 ,C-l 模 态 叶片 量 纲 归 一 化 振 型 的 
变形 程度 逐渐 减弱 ,C-2 模 态 的 叶片 变形 程度 也 逐 
渐 减 弱 , 但 变化 幅度 较 小 。 

图 10 给 出 了 为 1 ~5 情况 下 前 2 阶 耦 合 模 态 
频率 相对 偏差 的 变化 曲线 。 由 图 10 可 以 看 出 : 随 着 
叶片 数目 增加 ,前 2 阶 看 合 模 态 频 率 相对 偏差 逐渐 
降低 ; 随 着 半径 比 由 0. 10 增加 至 0.25,C-l 模 态 相 
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对 偏差 变化 曲线 依次 降低 ,而 C2 模 态 相对 偏差 降 
低 趋势 趋 于 平缓。 
1.0 


以 增加 


0 02 04 06 08 10 0 "02 G4 "08 GB ID 
叶片 轴 向 位 置 WP 时 片 轴 向 位 置 x/， 
(a) C-1 模 态 (b) C-2 模 态 


图 9 不 同 叶片 数 下 前 2 阶 稳 合 模 态 叶片 
量 纲 归 一 化 振 型 (c =0. 1 ,n=10’) 
Fig.9 ”Normalized first 2 coupling mode shapes of blade 


with various blade numbers (og =0.1 从 =10”) 
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(b) C-2 模 态 


(a) C-1 横 态 
图 10 不同 叶片 数 下 前 2 阶 耦 合 模 态 频率 相对 偏差 (人 =10 ) 


Fig. 10 Relative error of first 2 coupling mode frequencies 


with various blade numbers (人 =105) 

从 叶片 模 态 振 型 的 变化 角度 分 析 , 随 着 叶片 数 
目 增加 ,C-1 与 C2 模 态 的 振 型 弯曲 程度 均 有 所 降 
低 , 科 氏 效应 因此 减弱 ,相对 偏差 逐渐 减 小 ;同时 ,由 
于 半径 比 的 增加 降低 了 C-1 模 态 叶片 量 纲 归 一 化 振 
型 的 变形 程度 ,在 半径 比较 大 时 , 随 着 叶片 数目 增 
加 ,频率 相对 偏差 进一步 降低 。 与 之 相反 ,半径 比 增 
加 ,C-2 模 态 叶片 量 纲 归 一 化 振 型 变形 程度 增 大 , 与 
叶片 数目 对 振 型 的 影响 相互 抵消 ,因此 ,在 半径 比较 
大 时 ,C-2 模 态 频率 相对 偏差 基本 不 受 叶片 数目 影 
啊 。 上 述 结果 表明 ,多 叶片 - 浆 载 - 轴 耦 合 系统 受 科 
氏 效应 的 影响 相对 较 少 ,具有 更 为 良好 的 稳定 性 。 


3 结 论 


本 研究 利用 几何 精确 方法 结合 阶 次 截断 描述 时 


片 非 线性 变形 ,引入 了 科 氏 效应 的 影响 ,建立 了 叶 
片 - 桨 载 - 轴 耦 合 系统 动力 学 模型 ,根据 Hamilton 原 
理 推导 出 了 系统 振动 微分 方程 。 分 析 了 科 氏 效应 对 
耦合 模 态 的 作用 机 理 。 讨 论 了 不 同 叶 片 - 轴 刚 度 比 、 
叶片 - 桨 慌 羊 径 比 叶片 数目 下 科 氏 效应 对 系统 动力 
学 特性 的 影响 ,并 从 叶片 量 纲 归 一 化 振 型 的 角度 进 
行 了 机 理 分 析 。 

叶片 - 奖 载 - 轴 系 统 模 态 包含 叶片 弯曲 - 轴 扭 转 耦 
合 模 态 以 及 叶片 弯曲 非 耦合 模 态 。 科 氏 效 应 对 耦合 
模 态 存在 明显 影响 。 随 着 半径 比 的 增加 , 科 氏 效应 
对 C-l 模 态 频率 的 影响 逐渐 降低 , 而 对 C2 模 态 频 
率 的 影响 逐渐 升 高 ;刚度 比 在 10 ~ 10 时 ,发 生 了 叶 
片 模 态 转换 现象 , 模 态 转换 后 , 科 氏 效应 的 影响 显著 
增强 ;叶片 数目 增加 后 ,叶片 模 态 振 型 变形 程度 减 
小 科 氏 效应 的 影响 减弱 。 
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